132,962 research outputs found
cuIBM -- A GPU-accelerated Immersed Boundary Method
A projection-based immersed boundary method is dominated by sparse linear
algebra routines. Using the open-source Cusp library, we observe a speedup
(with respect to a single CPU core) which reflects the constraints of a
bandwidth-dominated problem on the GPU. Nevertheless, GPUs offer the capacity
to solve large problems on commodity hardware. This work includes validation
and a convergence study of the GPU-accelerated IBM, and various optimizations.Comment: Extended paper post-conference, presented at the 23rd International
Conference on Parallel Computational Fluid Dynamics (http://www.parcfd.org),
ParCFD 2011, Barcelona (unpublished
Removing the Stiffness of Elastic Force from the Immersed Boundary Method for the 2D Stokes Equations
The Immersed Boundary method has evolved into one of the most useful
computational methods in studying fluid structure interaction. On the other
hand, the Immersed Boundary method is also known to suffer from a severe
timestep stability restriction when using an explicit time discretization. In
this paper, we propose several efficient semi-implicit schemes to remove this
stiffness from the Immersed Boundary method for the two-dimensional Stokes
flow. First, we obtain a novel unconditionally stable semi-implicit
discretization for the immersed boundary problem. Using this unconditionally
stable discretization as a building block, we derive several efficient
semi-implicit schemes for the immersed boundary problem by applying the Small
Scale Decomposition to this unconditionally stable discretization. Our
stability analysis and extensive numerical experiments show that our
semi-implicit schemes offer much better stability property than the explicit
scheme. Unlike other implicit or semi-implicit schemes proposed in the
literature, our semi-implicit schemes can be solved explicitly in the spectral
space. Thus the computational cost of our semi-implicit schemes is comparable
to that of an explicit scheme, but with a much better stability property.Comment: 40 pages with 8 figure
A unified operator splitting approach for multi-scale fluid-particle coupling in the lattice Boltzmann method
A unified framework to derive discrete time-marching schemes for coupling of
immersed solid and elastic objects to the lattice Boltzmann method is
presented. Based on operator splitting for the discrete Boltzmann equation,
second-order time-accurate schemes for the immersed boundary method, viscous
force coupling and external boundary force are derived. Furthermore, a modified
formulation of the external boundary force is introduced that leads to a more
accurate no-slip boundary condition. The derivation also reveals that the
coupling methods can be cast into a unified form, and that the immersed
boundary method can be interpreted as the limit of force coupling for vanishing
particle mass. In practice, the ratio between fluid and particle mass
determines the strength of the force transfer in the coupling. The integration
schemes formally improve the accuracy of first-order algorithms that are
commonly employed when coupling immersed objects to a lattice Boltzmann fluid.
It is anticipated that they will also lead to superior long-time stability in
simulations of complex fluids with multiple scales
Hybrid finite difference/finite element immersed boundary method
The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian
description of the structural deformations, stresses, and forces along with an Eulerian description of the
momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary
methods described immersed elastic structures using systems of flexible fibers, and even now, most
immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This
work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian
variables that facilitates independent spatial discretizations for the structure and background grid. This
approach employs a finite element discretization of the structure while retaining a finite difference scheme
for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively
contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases
in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors
that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes.
The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse
structural meshes with the immersed boundary method. This work also contrasts two different weak forms
of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations
facilitated by our coupling approach
- …